Photon and neutral pion production in Au+Au collisions at = 130 GeV 
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We report the first inclusive plioton measurements about mid-rapidity (|j/j<0.5) from 
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Au+ Au collisions at ^/s^ = 130 GeV at RHIC. Photon pair conversions were reconstructed 
from electron and positron tracks measured with the Time Projection Chamber (TPC) of the STAR 
experiment. With this method, an energy resolution of AE/E ~ 2% at 0.5 GeV has been achieved. 
Reconstructed photons have also been used to measure the transverse momentum (pt) spectra of 
tt" mesons about mid-rapidity (|i/|<l) via the vr" — > 77 decay channel. The fractional contribution 
of the ■k'^ 77 decay to the inclusive photon spectrum decreases by 20% ± 5% between pt = 
1.65 GeV/c and pt — 2.4 GeV/c in the most central events, indicating that relative to tt'^ — > 77 
decay the contribution of other photon sources is substantially increasing. 



Relativistic heavy ion collisions provide the opportu- valuable probes of the dynamics and properties of the 
nity to excite matter into extreme conditions in the lab- resulting systems ^ JL Jsi i .fij ■ Unlike hadrons, which 
oratory. Of the particles which emerge from these col- have large interaction cross sections in dense matter, pho- 
lisions, photons are considered to be one of the most tons only interact electromagnetically and consequently 

have a long mean free path. This path length is typically 
much larger than the transverse size of the matter created 

in nuclear collisions j^. Therefore, with high probabil- 
* URL: |«wM . star ■ bni ■ gov| ity^ photons will escape from the system undisturbed, re- 
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taining information about the physical conditions under 
which they were created. 

Photons areproduced in all stages of heavy ion colli- 
sions 0, H S 0, II H ES in 111 Emi, from the first 
instant when the quarks and gluons of the opposing nuclei 
interact, through to long lived electromagnetic decays of 
final state hadrons. The production rate of photons dur- 
ing various stages of the created system has been theoret- 
ically calculated for a variety of initial conditions and sce- 
narios. It has been demonstrated that the emission rate 
from a hadron gas can be comparable to that expected 
from quark-gluon Compton and quark-antiquark annihi- 
lation processes in a net-baryon free system of deconfined 
quarks and gluons 7J. However, recent two-loop calcu- 
lations which include the quark bremsstrahlung process 
predict that photon production rates in quark matter ex- 
ceed those indicated by former one-loop calculations that 
only account for the Compton and annihilation processes 
[13. These calculations indicate "that the emissions 
from quark matter can outshine those from the hadronic 
matter" and that near and above pt = 2 GeV/c the 
contribution from hard scattered partons becomes more 
abundant than thermal photons from a hot hadronic gas 
at RHIC energy • Theoretical calculations and predic- 
tions like these underscore the importance of measuring 
photon spectra across a wide range of pt to investigate 
the matter created in heavy ion collisions. The measure- 
ments presented in this paper extend above and below 
the interesting region around = 2 GeV/c. 

Of all photon production mechanisms, late stage elec- 
tromagnetic decays of hadrons are the dominant source. 
At CERN SPS energies, photons from 7r° and r] decays ac- 
count for ^97% 5] of the inclusive photon spectrum. The 
remaining 3% arises from a combination of other sources, 
including electromagnetic decays of other hadrons such 
as the (jj, T]' and S". For these particles, thermal models 
that describe hadron production must be used to esti- 
mate the yields since their production rates have not yet 
been measured in heavy ion collisions at these energies. 
Measurement of the rj' p7 and S° Aj decays ap- 
pears to be promising with the energy resolution afforded 
by the photon reconstruction technique described in this 
paper. By estimating or measuring the yields of such par- 
ticles, their contribution to the single photon spectrum 
from electromagnetic decays can be calculated. However, 
the precision necessary to disentangle the rate of direct 
photon production from the rate of photons produced by 
electromagnetic decays is an experimental and theoreti- 
cal challenge. 

STAR has begun to address this challenge by mea- 
suring the spectra of both photons and tt^s. Pho- 
tons were measured by reconstructing pair conversions, 
"fZ—fe~e~^Z, with the electron and positron daughters 
detected in the STAR TPC. Reconstructed photons were 
used in turn to measure the rate of tt^ production via the 
— > decay channel. This paper discusses the tech- 
niques employed and the resulting spectra of photons and 
tt'^s that were measured. A full discussion of all cuts and 



variables used in this analysis may be found in 



I. DATA ANALYSIS 

The data presented in this paper were recorded by the 
STAR collaboration during the first ^/s^ =130 GeV 
Au-f Au run at RHIC. Events that had a primary colli- 
sion vertex position less than 100 cm and 150 cm distant 
from the geometric center of the TPC along the beam 
axis (z-axis) were selected for the photon and tt^ analyses 
respectively. Details of the STAR geometry are presented 
in Fig. n and discussed in 0, . For the event sample 
used for tt'^ measurements, which were limited by statisti- 
cal uncertainties, approximately 87% of the events in the 
minimum bias (least trigger-biased) data set passed the 
z vertex requirement. With this range of collision ver- 
tices, part of the support structure for the Silicon Vertex 
Tracker (SVT) could be utilized as a converter, as shown 
in Fig. n In addition, the inner field cage and gas (10% 
CH4 and 90% Ar) of the STAR TPC were also used as 
converters. The combined material from both detectors 
resulted in an average conversion probability of approx- 
imately 1% during the data run of year 2000. Although 
this conversion probability was low, it was compensated 
by the complete 2tt azimuthal acceptance of the STAR 
TPC. 

As discussed in j^, the definition of collision central- 
ity was based on the number of reconstructed primary 
tracks in the pseudo-rapidity range |77|<0.75. Using this 
as a basis, four centrality classes were defined, common 
to both the photon and 7r° analyses. They were an inclu- 
sive minimum bias (0-85% of the total inelastic hadronic 
cross section), peripheral (34%-85%), mid central (34%- 
11%) and central (0-11%). These centrality classes were 
selected to allow the extraction of tt" yield over a wide 
range of (0.25<pt<2.5 GeV/c) in independent regions 
of centrality. They contained 328980, 198196, 87484 and 
449095 events, respectively. 



A. Reconstructing Photon Pair Conversions 

The dominant interaction process for photons with 
a total energy above 10 MeV is pair conversion, 
'yZ^e^e^Z (Fig. |2J). Pair conversions that occurred 
in the detector material before or inside the TPC track- 
ing volume were reconstructed from the resulting charged 
particle daughters detected in the TPC. This was ac- 
complished in three steps: the selection of track, pair, 
and primary photon candidates. All three steps utilized 
the unique topological signature of a photon conversion - 
two tracks of opposite charge emerging from a secondary 
vertex with a small opening angle (w nieC^ /E-^ radians, 
where rUe is the electron mass and c is the speed of light). 

At the track level, improbable conversion daughters 
were removed by requiring tracks to satisfy a geometric 
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FIG. 1: Top figure: layout of the STAR experiment. Lower 
figures: density profiles of photon conversion points show the 
layout of the detector material. The structure at r = 46.5 cm 
is the inner field cage of the cylindrical TPC, while below r 
— 40 cm the SVT support cones and material are apparent. 



cut and to have the ionization energy loss expected for 
an electron in the TPC gas. Neglecting resolution ef- 
fects, the projection of a daughter track from a primary 
photon conversion onto the bend plane will form a circle 
which docs not enclose the collision vertex. This is be- 
cause photons typically propagate some distance before 
conversion and daughters emerge with a near zero open- 
ing angle (see Fig. Thus, low pt (<0.3 GeV/c) tracks 
with circular projections that enclosed the collision ver- 
tex in the bend plane of the 0.25 T solenoidal magnetic 
field were immediately removed. This cut was important, 
since the elimination of non electron (positron) tracks in 
this region of pt via ionization energy loss dE/dx is diffi- 
cult due to the fact that the highly populated pion band 
crosses the electron dE/dx band. It was not necessary 
to use this cut at higher pt, since the yield of particles 
drops and electron (positron) identification via dE/dx 
improves. It is also the case that at higher pt this cut be- 
gins to remove daughters of primary photons since stiff 
track geometries make the distance from the collision ver- 
tex to the closest point on the circular projection of the 
helix comparable to the resolution of the measurement. 
At all momenta, electron and positron candidates which 
had a dE/dx value between -2 and 4 standard deviations 
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FIG. 2: Schematic diagram of 7Z— s-e Z 



(ares) of the value expected for electrons and positrons 
were retained {ares denotes the resolution of charge par- 
ticle dE/dx measurements in the TPC gas; ares ~ 8.2% 
of the dE/dx value measured with a clean sample of elec- 
trons and positrons). The predicted energy loss curves 
for electrons, pions, kaons, and protons are shown as a 
function of rigidity in Fig. |21 The dE/dx requirement 
was chosen to be asymmetric, because on the lower side 
of the electron band other particle bands run in paral- 
lel and contamination is more prominent. On the upper 
side of the electron band other particle bands approach 
and cross the electron dE/dx band in a narrow range of 
momentum. This reduced the usefulness of a tight cut 
on the upper side. It was estimated that approximately 
3% of the true photon daughters were removed with this 
dE/dx requirement. 

Photon candidates were found by searching for track 
pairs which exhibited the topological signature of a pho- 
ton conversion. Oppositely charged tracks were paired 
and passed through a geometric filter. The filter required 
each pair to originate from a secondary vertex with a 
near zero opening angle and low invariant mass. Sec- 
ondary vertices were located by extrapolating daughter 
candidates to a common point. At the point of closest ap- 
proach, daughters were required to come within 1 cm of 
each other in the non-bend plane (rz-plane) and within 
1.5 cm of each other in the bend plane (xy-plane). The 
angular resolutions of opening angle measurements also 
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FIG. 3: Measured ionization energy loss in the TPC gas 
dE/dx versus rigidity. Highlighted curves indicate the pre- 
dicted dE/dx curves for electrons, pious, kaons, and protons. 
Photon selection criteria have not been applied to track can- 
didates in this figure. 
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FIG. 4: Invariant mass distributions of photon candidates 
assuming the electron and positron mass for the daughters. 
Left figure: the distribution can be separated into two con- 
tributions (indicated with dashed lines); a sharp lower-mass 
peak primarily composed of track-pair geometries that do 
not overlap in the bend plane, and a higher invariant mass 
peak composed of pair geometries that do overlap in the bend 
plane. Right figure: the invariant mass distribution of pho- 
tons calculated with the projection of the opening angle in 
the non-bend plane only. 



differed in the two planes. In the non-bend plane and 
bend plane the precision of opening angle measurements 
have single Gaussian sigmas near 0.02 and 0.1 radians, re- 
spectively. Since even at energies as low as 100 MeV pho- 
ton conversions on average have an opening angle of 0.01 
radians, ten times smaller than the precision in the bend 
plane, the full opening angle and the opening angle in 
the non-bend plane of each candidate pair were checked 
seperately. These values were required to be less than 
0.4 and 0.03 radians respectively. The differing angular 
resolutions in the bend and non-bend planes are appar- 
ent in the invariant mass distribution of pairs assuming 
an electron (positron) hypothesis for the daughters (Fig. 

The invariant mass distribution has a sharp peak 
near zero and a broad peak close to 0.012 GeV/c^. The 
sharp peak at lower mass primarily results from cases 
where the bend plane projection of the opening angle 
was assumed to be zero. In these cases, the track geome- 
tries do not overlap in the bend plane, so the tracks were 
assumed to be parallel in the bend plane at the point of 
closest approach. Track geometries that do overlap in the 
bend plane lead to a higher invariant mass because the 
complete opening angle was used in calculation of the in- 
variant mass. In this case, the less precise measure of the 
opening in the bend plane tends to dominate the result of 
the invariant mass calculation, moving and smearing the 
invariant mass peak. For this reason, a cut was placed on 
the value of the invariant mass of pairs calculated with 
only the non-bend plane projection of the opening angle. 
This cut required the invariant mass of candidate pairs 
to be less than 0.012 GcV/c^. The minimum mass re- 
turned by the calculation is 2me = 1.022 MeV/c^, which 
is above the first four 0.25 MeV/c^ wide mass bins in Fig. 
^ This causes the absence of entries in the lower mass 
bins in the figure. 

The kinematic parameters for photon candidates were 



derived from the kinematic variables of the associated 
daughter tracks. The energy was calculated by sum- 
ming the electron and positron energies. The angular 
direction in the bend plane was extracted by forming the 
cross-product of the vector from the helix center of the 
positron to the helix center of the electron, with the mag- 
netic field vector. The angular direction in the non-bend 
plane was found by averaging the direction of electron 
and positron at the conversion point. With these three 
variables all kinematic parameters for photon candidates 
could be derived. 

Primary photon candidates were selected from the set 
of all photon candidates by requiring the momentum vec- 
tor to be consistent with the direction of a photon orig- 
inating from the collision vertex. For these photons, the 
momentum vector has the same direction as the vector 
from the collision vertex to the conversion point. Due to 
differing angular resolution, the direction of the momen- 
tum vector and conversion point vector were compared 
separately in the bend and non-bend planes. For primary 
photon candidates, the difference between the momen- 
tum vector and the conversion point vector was required 
to be less than 0.035 and 0.015 radians in the bend and 
non-bend planes respectively. In order to reduce back- 
ground in the photon sample from the random pairing of 
primary tracks, conversion vertices in the region close to 
the collision vertex {rxy<lQ cm) were excluded. 



B. Photon Spectra 

Photon spectra were measured as a function of pt 
and y in three independent centrality classes as well as 
for an inclusive minimum bias data set. These spectra 
were produced from photon candidates identified with 
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the standard event, track, and photon selection criteria 
(discussed in Sec. 11 A|) . 

Photon yields were extracted using the particle identifi- 
cation information of the positive daughters in the TPC. 
A dE/dx deviant variable was constructed by compar- 
ing the energy loss predicted as a function of momentum 
for electrons and positrons with the measured dE/dx and 
rigidity of daughter candidates, folding in the dE/dx res- 
olution of the TPC. This variable accounts for the mo- 
mentum dependence in dE/dx, and its value is therefore 
independent of the daughter particle's momentum and 
the parent photon's pt- Consequently, the dE/dx de- 
viant values for daughters of differing momenta could be 
merged into bins based on the parent photon's pt and 
y. Distributions of the dE/dx deviant values of positive 
daughters were chosen rather than the negative daugh- 
ter to reduce the number of false photon candidates in 
the distributions that arise from knock-out electrons that 
originate when charged particles scatter in the detec- 
tor material (i5-electrons). The remaining contamination 
from this scattering process, which may result in a knock- 
out electron and positive particle (tt, K, or p) having a 
momentum in a region where the dE/dx bands overlap 
the positron band (see Fig. |3Jl, was removed by requiring 
the fraction of the positive daughters energy to the total 
photon energy to be less than 75%. The shape of the 
remaining background in the dE/dx deviant distribution 
was studied on a sample of photon candidates that satis- 
fied anti-photon cuts (primarily a sample of background 
candidates). Anti-photon cuts suppress positrons from 
true photons by requiring the photon selection criteria 
to be in the outermost extent of the cut distributions. 
For example, the two-track distance of closest approach 
for daughter tracks was required to be 1.5<|da;i;|<2 cm 
and l<|dz|<1.5 cm. A two parameter exponential plus 
linear function was used to describe these background 
distributions, as shown for one pt bin in Fig. [S] Param- 
eters of the background functions were found by fitting 
to the background distributions with the region around 
the expected value removed {dE/dx deviants between - 
1.5 and 3 ares)- This was necessary to avoid fitting the 
signal from residual photons that still existed after the 
application of anti-photon cuts. With the knowledge of 
the background shape, the raw yield of photon candidates 
was extracted using a three parameter Gaussian function 
plus the background function which had one free scaling 
parameter (also shown in Fig. [SJ. 

The purity of the photon candidate sample was deter- 
mined by dividing the integral of the Gaussian function 
by the integral of the entire Gaussian plus background 
function between dE/dx deviant values of -2 and 4 ares- 
For Pt < 0.75 GeV/c, the purity of the photon candidate 
sample is greater than 90% in all centrality classes. In the 
0-11% most central centrality class, where the purity is 
the lowest, the purity drops linearly from approximately 
90% at Pt = 0.75 GeV/c to about 60% at pt = 2.4 GeV/c. 
A cleaner sample (purity >95% below pt = 0.90 GeV/c 
for the 0-11% most central collisions) was obtained by 



requiring photons to convert in the inner field cage and 
TPC gas, ra;y>40 cm. 

Uncorrected yields were obtained from the weighted 
sum of the entries in the distributions. The weights were 
extracted by dividing the height of the Gaussian func- 
tion by the height of the entire fit at the location of each 
entry. Distributions in y were extracted in a 4 x 10 ar- 
ray of Pt~y bins to properly account for the variation in 
efficiency as a function of pt- Three 0.25 GeV/c wide 
Pt bins were used below pt — 0.75 GeV/c where the ef- 
ficiency grows rapidly, and one large bin was used for 
0.75<pt<2.5 GeV/c where the efficiency is flat. The pt 
distributions did not require division into pt and y bins, 
since the corrected y distributions and the input distri- 
butions for the efficiency calculations were both uniform 
in y for |?;|<0.5. 



Background 0.9<P,<1 .05 GbV/C yield = 25579 +/- 198 




dE/dx deviant ( a) dE/dx deviant ( a) 



FIG. 5: dE/dx deviant distributions of positive daugh- 
ters with 0.9<pt<1.05 GeV/c from photon candidates with 
rxy>lQ cm. Left figure: distribution of background candi- 
dates fit with an exponential plus linear function. Right fig- 
ure: distribution of the positron signal fit with a Gaussian 
function plus a scaled background function. 

Efficiency corrections were applied to each pt-y bin 
independently. These corrections were calculated with 
detailed simulations (GEANT 3.21) of the propagation 
of photons and daughter particles through a realistic de- 
tector geometry. A TPC Response Simulator (TRS) was 
used to simulate the drift and electronic response of ion- 
ization deposited in the TPC. Digital pad signals pro- 
duced by TRS were embedded pixel-by-pixel into real 
events. 

Each simulated event contained 2000 photons gener- 
ated flat in Pt and y. On average, approximately 20 
of these photons interacted with the detector material 
(a consequence of the low conversion probability). This 
added the ionization of about 40 daughter tracks to each 
real event. This number is less than 2% of the number 
of charged particles in the embedded phase space in a 
typical high multiplicity event. Therefore even in high 
multiplicity events, which are most sensitive to over- 
embedding, the introduction of 2000 photons into each 
event had a negligible effect on the track reconstruction 
efficiency. An association process was used to link re- 
constructed and generated photons. The photon finding 
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efficiencies for different centrality definitions were calcu- 
lated by dividing the distributions of reconstructed pho- 
tons correctly associated with a generated photon by the 
input distributions of generated photons. 
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FIG. 6: Number of reconstructed photon conversions as a 
function of the conversion location for both real data (solid 
Une) and GEANT simulated events (dashed line). Left figure: 
conversion density as a function of radial distance from the 
beam axis. Right figure: conversion density as a function of 
the distance along the beam axis, z. 



To reveal systematic trends caused by differences in the 
layout of the detector material and the material map used 
in simulation to calculate efficiency corrections, spectra 
were produced with different requirements on the mini- 
mum distance between the location of conversion and the 
beam axis in the xy-plane [rxy>^^ cm and r^y>\Q cm). 
It was found that an additional correction factor was 
needed to compensate for differences in the layouts of 
the detector material between Vxy = 10 cm and r^y = 

= 40 cm, 



40 cm, as illustrated in Fig. El Above 
the material of the inner field cage and the gas of the 
TPC were well described in the simulation. This region 
was used as a reference to calculate the correction fac- 
tor (1.42) for the rxy>V^ cm spectra. All -pt and y data 
points in the rxy>y^i cm spectra were linearly scaled by 
this factor. 

Corrected pt and y spectra for the various centrality 
classes are shown in Fig. [7] Systematic uncertainties 
of 7% point-to-point and 12% overall (correlated) have 
been estimated for the measurements in the -pt spectra. 
There is a 12% (correlated) systematic uncertainty in the 
normalization of the y spectra. These uncertainties ac- 
count for uncertainty in the detection efficiency and po- 
tential measurement bias that may arise from differences 
in the physical and simulated material maps. 

Values in the lowest rxy>Wi cm pt bins, 
0<pt<0.15 GeV/c, are systematically 15-25% lower than 
in the corresponding ra;y>40 cm bins. This is attributed 
to the efficiency correction being underestimated in these 
bins for the rxy>\^ cm, because the simulation lacked 
material between 10<ra;y<40 cm. This gave the inner 
field cage and TPC gas a larger fraction of the total 
conversion probability in the simulation, and resulted 
in a mean r^y conversion point that is closer to the 
TPC. The combination of the shift in the mean position 
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Corrected photon pt (left) and y (right) spectra for 
~ = 130 GeV. The pt distribu- 
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Au-|- Au collisions at 



tions are for mid-rapidity photons, |2/|<0.5. Statistical uncer- 
tainties are shown. Systematic uncertainties in the pt spectra 
and on the normalization of the dN/dy spectra have not been 
included in this figure. 



of conversions and the linear scaling of the spectra 
to compensate for differences in the material layouts 
artificially increased the efficiency of low pt photons in 
the rxy>\Q cm spectra. Therefore the corrections, which 
are the inverse of the efficiencies, are too small in the 
ra;j,>10 cm spectra at low pt. 



C. Measuring the yield of tt" — > 77 decays 

The uncorrected yield of tfi mesons was extracted from 
the invariant mass distributions of photon pairs in var- 
ious "Pt bins. Individual decays could not be uniquely 
identified, because of the large combinatorial background 
in the invariant mass distributions. The combinatorial 
backgrounds were simulated by rotating the momentum 
vector of one photon in each pair by tt radians in the 
bend plane. In this way, it was possible to create combi- 
natorial background distributions which preserved event 
characteristics such as the vertex position along the beam 
axis, the event multiplicity and anisotropic flow. Due to 
the azimuthal symmetry of the STAR TPC, this type 
of rotation also ensured that a consistent geometric ac- 
ceptance and track reconstruction efficiency were main- 
tained. At the same time these rotations moved and 
smeared the invariant mass values of the pairs that are 
correlated through two photon decays. The shape of 
the resulting background distributions near the 7r° mass 
(±0.1 GeV/c^) was well described and smoothed with a 
second order polynomial (see Fig. |SJ). This functional 
form was also used to describe the shape of the combina- 
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torial background in the unrotated invariant mass distri- 
butions. A Gaussian function was used to describe the 
enhancement at the tt" mass. The width (cr) of the en- 
hancement ranged between 4 and 15 MeV/c^, and was 
found to be consistent with simulations in all central- 
ity classes and function of pt- It was dominated 
by the photon momentum resolution not the intrinsic 
mass width of the 7r° (~8 eV/c^). For this reason a 
Breit-Wigner function that would be appropriate to de- 
scribe the intrinsic width of such a resonance was not 
used. Equation ^ is the complete function that was 
used to describe invariant mass distributions of photon 
pairs. Systematic uncertainty related to the choice of the 
invariant mass bin width, 5, was studied by comparing 
the measured yields for various choices of 5. These stud- 
ies indicated that the systematic uncertainty due to this 
effect was much smaller than statistical uncertainties in 
the measurements. 



C {x=M^^) = +B(a + bx + cx^) , 

(TV27r 

where S is the width of the invariant mass bins, 
N is the number in the Gaussian peak, and B is 
the scale factor of the background function. 

A feature of this method of photon reconstruction is 
that the location of the enhancement from tt" decays 
in the two-photon invariant mass distribution is a few 
MeV/c^ lower than expected (see Fig. IS)). This is at- 
tributed to energy loss experienced by the electrons and 
positrons in the detector material. The "global" tracking 
routine used in this analysis only compensated for energy 
loss in the TPC gas and not for that in other detector 
material. This resulted in the reconstructed momentum 
for electrons and positrons originating prior to the gas 
volume of the TPC to be systematically lower than their 
original momentum. The small ('--^l MeV on average) en- 
ergy loss experienced by each of the 4 daughter particles 
translated to a few MeV/c^ shift in the location of the tt" 
invariant mass peak. This hypothesis is consistent with 
a similar feature in simulated events. The location of the 
reconstructed invariant mass peak for simulated n^s sys- 
tematically decreased as the radial distance between the 
beam axis and conversion point of the closer photon in 
the pair decreased. This implies that a larger 7r° mass 
deviation occurs when more detector material is between 
the conversion point and the TPC. The reconstructed en- 
ergy of simulated photons was also systematically lower 
than the energy input to the simulation and larger de- 
viations also occurred when more detector material was 
between the conversion point and the TPC. 



D. Spectra of tt mesons 

Two iterations (described in Sec. II C|l were performed 
to fit the data and extract the yield of tt^s as a function 
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FIG. 8: Two-photon invariant mass distributions for can- 
didates with 0.75<pt<l GeV/c in the 0-11% most central 

197 197 

Au-I- Au collisions at ^/s^ = 130 GeV. Top frame: in- 
variant mass distribution with one photon rotated by tt ra- 
dians fit with a second order polynomial. Middle frame: in- 
variant mass distribution of photon pairs fit with Eq. Q; 
the background function is also shown. Bottom frame: in- 
variant mass distribution after the combinatorial background 
was removed. The enhancement near the tt" mass is located 
at 0.131 GeV/c^ and has a Gaussian sigma of 0.006 GeV/c^ 



of Pt- For the first iteration, 4 free parameters were used 
in the fit: the yield (N), mass (m), and width (cr) of the 
Gaussian function describing the peak, and the scale fac- 
tor for the background function (B). For the second iter- 
ation it was assumed that the width of the invariant mass 
peak increased linearly with pt, as seen in the simulated 
events. The width parameters for the fits in the second 
iteration were obtained from a linear fit to the widths 
found in the first iteration. The values of the slopes for 
these linear fits (« S^AIeV/c^) / {GeV/c)) were consis- 
tent with those found in simulation. The width param- 
eters were then fixed to the value of the linear function 
at the center of each pt bin. This reduced the number of 
free parameters in the second pass to three and increased 
the stability of the fits. 

Uncorrected yields about mid-rapidity (|?;|<1) were 
extracted in various pt bins for the four different cen- 
trality classes. The narrow 6 MeV/c^ width (sigma) of 
the enhancement at the tt'-' mass measured with this re- 
construction method is significantly better than what is 
typically achieved using a conventional lead-scintillator 
sampling calorimeter (20 MeV/c^ sigma). The improve- 
ment is a result of the excellent photon energy resolution 
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(3% at 1 GeV) obtained with this method of photon re- 
construction. The narrow width improves the signal to 
background ratio and enables the extraction of raw 7r° 
yields at low pt (pt<0.75 GeV/c) where the signal to 
background ratio is seriously degraded by a large combi- 
natorial background. 

Efficiency corrections were calculated with a procedure 
similar to the one used to calculate the photon detec- 
tion efficiencies (described in Sec. IIB|I . except that 7r°s 
were selected in GEANT. Only the ionization in the TPC 
gas from daughters of those tt^s selected in GEANT was 
passed to TRS. This was necessary to perform the cal- 
culation in a reasonable amount of cpu time and to pre- 
vent saturating real events with ionization. One conse- 
quence of the low conversion probability is that on av- 
erage only 1 in 10000 tt" 77 decays is expected to 
be detected through the reconstruction of pair conver- 
sions. At the same time about 1 in 50 7r° 77 de- 
cays produces a photon that converts to create a pair 
of tracks that ionize the gas in the TPG. Selecting de- 
tectable tt" — > 77 decays in GEANT (tt^s that decay into 
two photons, with both photons undergoing interactions 
which create at least one daughter within the TPC accep- 
tance) reduced the amount of uninteresting ionization in 
each simulated event. With this selection, the insertion 
of up to 12 detectable tt" 77 decays into each event 
added less than 2% to the number of tracks in the phase 
space of the embedding. These events were reconstructed 
with the same software used to reconstruct real events. 
Reconstructed photons that could be associated with a 
photon from a simulated tt^ — > 77 decay were retained. 
These photons were used to generate two-photon invari- 
ant mass distributions. The yields of simulated tt'^s were 
calculated in the same pt and centrality bins, and with 
the same fitting procedure as the raw yields. Efficiency 
corrections for these bins were obtained by dividing the 
reconstructed distributions by the input distributions. 

Corrected pt spectra of tt^s were obtained by applying 
efficiency corrections to the pt distributions of the raw 
yields. The corrected spectra are shown in Fig. |^ The 
uncertainties shown are statistical and mainly reflect the 
low number of real tt^s measured. They combine the 
uncertainty in the raw yields and efficiency corrections. 
Systematic uncertainties due to the cuts used were stud- 
ied by varying track, photon, and tt*' cuts. These stud- 
ies revealed that the statistical fluctuations dominate the 
systematic variations in this analysis. 

Corrections to the normalization of the spectra were 
made to compensate for the overall difference in the pho- 
ton conversion probability in the real detector material 
and that used in simulation. These correction factors 
were based on the corrected yield of 7r°s (Yr^j^>4o) for 
photons that converted in the inner field cage or the TPC 
gas {rxy>4:0 cm) where the two material maps are con- 
sistent with each other (see Fig. (H)). For l<pt<2 GeV/c, 
a centrality-independent material correction factor was 
directly extracted by dividing Y™J1^4q by Y™'" for min- 
imum bias events. In this pt interval, the ratio of the 
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FIG. 9: Transverse momentum spectra of the tt" meson about 
mid-rapidity, |j/|<l, for different centrality bins. Along with 
the point-to-point statistical uncertainties shown, the spectra 
have a common uncertainty in the normalization of ±40%. 
The central data (0-11%) has an additional normalization 
uncertainty of ±19%, which results in a total normalization 
uncertainty of ±49%. Dashed lines indicate Boltzmann fits 
to the spectra. 



tt" efficiency for Y™'^4q to the efficiency for Y"*™ is uni- 
form. Therefore it was not necessary to correct this factor 
for interplay between the shape of the efficiency correc- 
tion and the exponentially falling spectrum. This mate- 
rial correction factor was crosschecked in three different 
Pt windows. The resulting variation in the the factors 
was found to be 5 x smaller than the uncertainty in an 
individual factor (~40%). This uncertainty is statistical 
and arises due to the relatively small number (~50) of 
tt'^s reconstructed from photons with r2;y>40 cm. The 
stability of this measurement was checked by compar- 
ing the results for different invariant mass bin widths. 
This check confirmed that the statistical uncertainty in 
the raw yields is the dominant uncertainty. The mate- 
rial correction factor for the minimum bias data set was 
0.31 ± 0.12, with l<pt<2 GeV/c as the window of 7r° pt. 
This factor was used to scale the normalization of cen- 
trality classes formed by taking subsets of the minimum 
bias triggered data set (0-85%, 34%-85%, 34%-ll%). 
The uncertainty in this factor (±0.12 or 40%) is common 
for all these centrality bins and cancels out when ratios 
are taken. This factor can not be used to normalize the 
0-11% centrality bin, because of the differing z vertex 
distributions between the central and minimum bias trig- 
gered data sets. A separate centrality-independent fac- 
tor was calculated for \z^ertex\<'^^ cm where the vertex 
distributions are similar in the two data sets. Minimum 
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bias triggered events with z vertices in this region were 
used to calculate the factor. The factor, 0.26 ± 0.11, 
was obtained by dividing Y™'!>4q by Y™"^|^.^g. Using 
this factor, the normalized yield for the 0-11% central- 
ity bin was computed by taking the product of the fac- 
tor and Y|=^^^"^*^|^75. This normahzed yield, Y^'^"* = 12.0 
± 5.7 for l<pt<2 GeV/c, is independent of the vertex 
distribution. The final 0-11% spectrum was scaled by 
Ycent/Y|=|^"j l^^gg = 0.27 ± 0.14. In summary, the 0-11% 
most central class of events has the same common uncer- 
tainty of ±40% plus an additional uncertainty of ±19%. 
These two uncertainties were combined (±49%) for the 
purpose of comparison to other spectra, like the spectrum 
of charged hadrons. 

II. CONTRIBUTION OF THE tt" ^ 77 DECAY 
TO THE INCLUSIVE PHOTON SPECTRUM 

Electromagnetic decays of neutral mesons are the dom- 
inant source photons in heavy ion collisions. Among 
these, the 7r° 77 decay is largest contributor to the 
spectrum of inclusive photons. Its large contribution 
hides the signal from other sources, such as direct pho- 
tons emitted during the early stages of heavy ion colli- 
sions. 

To investigate how the tt'^ — > 77 decay contributes to 
the inclusive photon spectrum, the distributions of tt^s 
were used to generate the single photon spectrum ex- 
pected for the daughters of 7r° 77 decays. The pt 
distributions of tt*^ were fit with both a Boltzmann func- 
tion and a Bose-Einstein function. For both functions, 
the total energy of the tt*^ was replaced by its transverse 
energy (^/pf ± m^c'*) under the assumption that the 
system is boost invariant near mid-rapidity. This as- 
sumption is supported by the flat shape of particle ra- 
pidity distributions close to mid-rapidity p3.l23j|. Other 
more sophisticated functions, incorporating resonances 
that decay into tt'^s and/or handling radial expansion of 
the system in more detail with additional parameters, 
were not chosen because the additional parameters were 
not well constrained by the 7 or 8 data points of the 
spectra. Both the Boltzmann and Bose-Einstein func- 
tions treat the system as a thermalized gas and converge 
to exponential functions at high pt ■ 

Distributions of the "pt dependence of tt'^s were gen- 
erated using these functions, assuming that the rapid- 
ity and azimuthal distributions are flat. The input ra- 
pidity distribution of tt^s was limited to |j/|<2. This 
rapidity window produces more than 99% of the pho- 
tons with |j;|<0.5 from vr" — > 77 decays. These distribu- 
tions were passed through a Monte Carlo decay simula- 
tor used to calculate the tt" — > 77 decay kinematics and 
boost between the center of momentum and laboratory 
frames. The momentum information of the decay pho- 
tons was used to produce the single photon -pt spectra of 
the daughters. 

The fraction of the photons from 7r° — > 77 decays in the 



inclusive photon spectra was calculated as a function of 
Pt by dividing the simulated spectra by the measured in- 
clusive photon spectra, as shown in FigureEB The shape 
of the resulting distributions for various centrality classes 
is independent of the uncertainty in the normalizations 
of the TT*^ spectra, but may depend on the assumed pt 
dependence of the tt" spectra. For this reason, the frac- 
tions are shown for pt>0.45 GeV/c, where the photon 
contribution is determined from tt'^s in and above the 
measured pt interval. The kinematics of tt*^ 77 decay 
limit the pt of the daughter photons. For example a ■n^ 
of pt<0.435 GeV/c can only produce photons with pt < 
0.45 GeV/c. Thus, the unmeasured portion of the tt" 
spectra below pt — 0.25 GeV/c does not contribute to 
the photon spectra in the region where the fractions are 
plotted. 

The fraction of photons from vr" — s- 77 decay in the 
inclusive photon spectrum is approximately constant be- 
tween 0.75<pi<1.65 GeV/c. For the 0-11% most cen- 
tral event class, the fraction begins to decrease substan- 
tially near pt = 1.65 GeV/c assuming either the Boltz- 
mann function or Bose-Einstein function to describe the 
Pt spectra of tt^s. Speciflcally, the relative contribution 
from 7r° 77 decays decreases by 20% ± 5% from pt = 
1.65 GeV/c to pt = 2.4 GeV/c. Both point-to-point un- 
certainties in the photon spectrum as well as substantial 
uncertainty in the slope parameters of the Boltzmann 
(0.281 GeV ± 0.013) and Bose-Einstein (0.289 GeV ± 
0.014) flts have been included in the 5% uncertainty in 
the region of the decrease. A similar trend was observed 

208 208 

by the WA98 collaboration in Pb+ Pb collisions at 
= 17.2 GeV HQ. The WA98 coUaboration re- 
ported an excess of photons above 1.5 GeV/c in cen- 
tral collisions after accounting for photons from all ex- 
pected electromagnetic decays. For electromagnetic de- 
cays other than the tt*^ —> 77 decay, the 77 — > 77 decay 
channel is expected to be the next largest contributor. 
From Pt = 1 GeV/c to pt = 4 GeV/c, its contribution 
is expected to be approximately 15% and to be fairly 
uniform in pt, increasing less than 5% per GeV/c. The 
WA98 collaboration has estimated that the summed con- 
tribution of all other electromagnetic decays is less than 
a few percent at SPS energies g]. Based on the above as- 
sumptions (the tt'^ Pt spectrum has a Boltzmann or Bose- 
Einstein Pt dependence, the 77 contribution increases by 
less than 5% per GeV/c, and the summed contribution 
of all other electromagnetic decays is less than a few per- 
cent) it is unlikely that electromagnetic decays fully ac- 
count for the observed single photon yields in the 0-11% 
most central event class. 



III. COMPARISONS TO PUBLISHED DATA 

Comparisons between the 0-11% 7r° spectrum and the 
0-10% charged hadron spectrum were used to study the 
composition of the hadron spectrum as a function of p(. 
The ratio of the 7r° data points to a power law function flt 
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FIG. 10: The Ratio as a function of pt of the distributions 
of photons from vr" — > 77 decays to the measured photon 
spectra. The tt" 77 photon distributions were generated 
assuming that the pt dependence of the tt" pt spectra follow 
a Boltzmann distribution. These ratios include both statisti- 
cal and systematic uncertainties in the photon spectra. Un- 
certainties in the normalization of the tt'^ spectra arise as pt 
independent uncertainties and have not been included. Nor- 
malization uncertainties in the 7r° spectra are 40%, correlated 
between all ratios, with an additional 19% uncorrelated un- 
certainty for the 0-11% centrality ratio. Uncorrelated pt de- 
pendent uncertainties that arise from the uncertainty in the 
slope parameters of fits to the tt" pt spectra (11%, 9% and 
5% respectively for the 34-85%, 11-34% and 0-11% central- 
ity classes) have not been included. 
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FIG. 11; Top frame: comparison between STAR vr" and in- 
clusive charged hadron 4- /i^) /2) spectrum f2f3 about 
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mid-rapidity for Au+ Au collisions at ^fs^ — 130 GeV. 
Bottom frame: ratios of these spectra to a power law fit to 
the STAR inclusive charged hadron spectra. For reference. 



dashed lines indicate 



of 0.5 ± 0.25 and 



1. Normalization uncertainties in both the tt'^ (±49%) and in- 
clusive charged hadron (±11%) measurements have not been 
included with the data points. 



to the charged hadron spectrum ((h~ + /2) is shown 
in Fig. ^2 At Pt = 2 GeV/c the ratio of 7r°s to charged 
hadrons approaches 50% (also shown in Fig. Ill|l . As- 
suming isospin symmetry for charged and neutral pions 
{{dN^+ + dN^-) = 2 (dN^o)) the proton to pion ratio is 
close to 1 at = 2 GeV/c. This result is similar to the 
previous observation in central collisions that the "p and 
p yields are comparable to the 7r+/~ yields." |2^. An- 
other method of probing the ratio of baryons to mesons 
in the system is by examining the ratio of tt" to A pro- 
duction (shown in Fig. I12|) . In this figure, Bose-Einstein 
functions have been used to describe the A and A spectra. 
The value of the 7r° data points has been divided by the 
value of the Bose-Einstein function describing the A and 
A spectra at the center of the bins to obtain the tt" to A, 
and TT^ to A ratios. At pt = 2 GeV/c these ratios are 
approximately 1, consistent with other measurements of 
the baryon to meson ratio. 

A comparison of the tt" spectrum for the 0-11% cen- 
trality class was also made to other identified pion spec- 
tra for central collisions at ^s„„ = 130 GeV. The 
PHENIX experiment has published 7r°, tt"*" and 7r~ spec- 



tra for central events 0, . The 7r° spectra were mea- 
sured via the n'^ —> 77 decay channel using both Lead- 
Scintillator (PbSc) and Lead-Glass (PbGl) calorimeters. 
These data overlap the STAR tt" measurement in the 
range l<pt<3 GeV/c. Ratios between the central STAR 
tt" spectrum to power-law fits of the PHENIX tt" spec- 
tra indicate that the shapes of the spectra are consistent 
(Fig. I13|l although the two experiments have a system- 
atic offset in normalization. In the region of overlap, 
the STAR 7r° spectrum is systematically higher than the 
PHENIX spectra. A systematic difference in the same 
direction is also observed in comparisons between the 
charged hadron spectra ((/i+ ± /2) from the two ex- 
periments for l<pt<3 GeV/c. Direct comparison of the 
0-11% STAR TT° spectrum and the 0-5% PHENIX 7r± 
spectrum shows the two are consistent in shape (Fig. 
I14|l . although once again the normalizations are system- 
atically different once a linear scale factor (0.91 ± 0.04, 
deduced from C values given in 25]) is applied to con- 
vert the 0-5% TT^ data to the 0-10% centrality class. 
These ratios indicate that meson to baryon ratios are in- 
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FIG. 12: Top frame: comparison between STAR tt", and A(A) 
measurements about mid-rapidity for Au+ An colli- 
sions at y^s^ = 130 GeV. Lower frames: ratios of these spec- 
tra to Bose-Einstein function fits to the STAR A and A mea- 
surements. For reference, dashed lines indicate dN / dN iyf^x)fit 
of 1. Normalization uncertainties in both the tt" (±49%) and 
A(A) (±10%) measurements have not been included with the 
data points. 



FIG. 13: Top frame: comparison between the STAR tt mea- 
surement and PHENIX tt" measurements 12^ about mid- 
197 197 

rapidity for Au± Au collisions at y^s^ = 130 GeV. 
Lower frames: ratios of these spectra to a power law fit to 
the PHENIX tt" measurements. For reference, dashed lines 
in the lower frames indicate dN/dNpg^^jx nOfit of li ''■nd 
1.78 ± 0.98 and 1.64 ± 0.86 in the middle and bottom frames 
respectively. Normalization uncertainties in both the STAR 
(±49%) and PHENIX (±25% for PbGl and ±20% for PbSc) 
tt" measurements have not been included with the data points. 



ternally consistent within PHENIX and STAR, although 
between PHENIX and STAR the normalization of the 
spectra is systematically shifted for pt between 1 GeV/c 
and 3 GeV/c. 

IV. CONCLUSION 

We have presented the first inclusive mid-rapidity, 

|y|<0.5, photon spectra as a function of centrality from 

197 197 

Au± Au collisions at = 130 GeV. The spec- 

tra of TT^s about mid-rapidity, |?;|<1.0, have been pre- 
sented; as well as the contribution from tt^ 77 decays 
to the inclusive photon spectrum. Near pt = 1.65 GeV/c 
the fractional contribution from tt'^ 77 decays to the 
inclusive photon spectrum for the 0-11% most central 
collisions begins to decrease significantly. This decrease 
indicates that relative to the 7r° 77 decay, the con- 
tribution from other sources of photons increases with 
Pt. In order to understand the origin of this decrease 
other electromagnetic decays must be measured or esti- 



mated. The combination of increased event statistics in 
future measurements with the excellent energy resolution 
achieved using this photon detection technique (AE/E w 
2% at 0.5 GeV/c) will make the measurement of the 77 
feasible. A statistically significant enhancement in the 
two-photon invariant mass distribution has already been 
observed in the vicinity of the 77 mass. Increased event 
statistics will also lead to higher precision measurements 
and extend the pt range of the 7r° spectra. Advances in 
these directions will not only enhance our understanding 
of contributions to the single photon spectra, but will also 
aid measurements of the relative abundance of mesons 
and baryons at high pt (>3 GeV/c) where the expected 
effects of collective motion become less dominant. 
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the tt" (±49%) and charged pion (±11%) measurements have 
not been included with the data points. 
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